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q ■ Abstract 

The neutron structure function F%{x) is evaluated within the kinematic range 
10~ 3 < x < 1 from the deuteron and proton data by employing relativistic the- 
H-j ! oretical description of F^^x) and several assumptions on the high-x asymptotics 

' of F^{x) / F^x). It is shown that new measurements of F^{x) in the range 0.6 < 

x < 0.8 would substantially improve understanding of the relation between d and u 
, valence quarks in the limit x — * 1. 

> : 

: 1 Introduction 

The valence quark structure of the proton and neutron has for some time been assumed 
to be understood. A number of experiments have provided a detailed representation of 
the nucleon's quark distributions over a wide range of kinematics with some exceptions: 
the range of Bjorken x close to a single nucleon kinematic limit, x = 1 , remains to be 
inaccessible. Actually, the situation is even more unfortunate because the range beyond 
x = 0.7 has been poorly explored experimentally. Recent reviews of the neutron and 
proton structure, which still do not pretend to be complete, are presented in Refs. [1, 2]. 
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The valence u and d quark distributions are generally obtained from measurements of 
the proton and neutron structure functions, F£ (x) and F^x), respectively. The u quark 
distribution is relatively well constrained by the F^x) data for x < 0.8 but the absence 
of free neutron targets has resulted in large uncertainties in the d quark distribution 
beyond x > 0.6. Major uncertainties arise from model considerations of the deuteron 
structure function F 2 D (x), from which Fg(x) is evaluated. They are normally represented 
as a spectrum of different ratios i? n//p = F^jF^ as a function of x. There is no theoretical 
constraint on the x dependence of i? n//p . Numerous very different constraints on R n / p 
at the kinematic boundary x — 1 have been suggested in quark models [3, 4] and QCD 
inspired models [5, 6]. They rely on model considerations of the d/u ratio in the limit 
x — > 1 and neglect the contribution from sea quarks. Such predictions offer a nice testing 
ground for our understanding of the role which valence quarks play in the nucleon wave 
function. However, it is extremely difficult to confront the predictions made for the limit 
x — *■ 1 with measured values which have to be extrapolated to this limit. This explains 
coexistence of numerous models of F®(x) used in practice of measurements of F£{x) as 
well as for motivating new experimental research into the valence quark structure. 
As it is demonstrated in Ref. [7] by using examples of evaluating the d/u ratio, im- 
provement of the knowledge of F^x) in the region x > 0.75 is very important for many 
applications in hadron physics. The procedure of extraction of F£(x) from the F®(x) data 
and its theoretical justification was a subject of many publications (see e.g. Refs. 



which defended a rigorous consideration of the role of nuclear effects. From the intense 
discussion in Ref. [9] one learns that compromises and simplifications distort consider- 
ably -F 2 n (x) and, therefore, d/u in the high x range. Theoretical understanding of nuclear 
effects is further discussed in Ref. [10] by presenting model independent relations which 
follow from the concept of quark-hadron duality. 

The present paper continues the series of publications [7] - [12] by suggesting an alternative 
approach to extraction of the F^x) from the data collected in deep inelastic scattering 
(DIS) experiments which relies on relativistic theoretical description of FJ?(x) [13, 14] and 
well defined assumptions on the high-x asymptotics for R n ^ p . We forward our criticism 
against unjustified simplifications frequently made in consideration of the "nuclear" effects 
in F 2 D . Most common misapprehension shows up in attempts to find an analogy (and even 
extrapolation rule) between the EMC effect and a modification of the nucleon structure 
inside the deuteron. There is no such theoretical concept as the EMC effect, which is 
just a bare observation that F^jF^ < 1 in a certain range of x. Therefore, there are no 
grounds to relate directly the difference between F 2 A and F^ with that of _F 2 D and the free 
nucleon structure function F 2 N , where F 2 N = (F 2 n + .F 2 p )/2. The difference between F 2 N and 
F 2 D can be conceptually original [15]. 

Experimental information on the deuteron and proton structure functions is available 
from the experiments of BCDMS, EMC, SLAC, E665, NMC, HI [16]. The data for F 2 D 
and F% and their ratio are available in the range 10~ 3 < x < 0.6. There is also data at 
0.6 < x < 0.9 from SLAC with relatively large errors. The SMC collaboration proposed 
overall fit of world data which gives F 2 and F 2 D in the range 10~ 4 < x < 0.85 [17]. We 
make use of this approximation and extend it up to x — 1 based on the quark counting rule 
for the x — > 1 limit of F 2 P . Further we consider extraction of Ff(x) by assuming that the 
deuteron can be considered as the two-nucleon bound state. We show that the behavior of 
F£(x) outside the x range covered by measurements can be established by employing well 
known theoretical prescriptions. Possible ambiguities connected with these prescriptions 
are investigated. 

2 Neutron structure function 

Our consideration of the deuteron is based on the standard picture of the proton and 
neutron bound together into the simplest nucleus. It is then a favored source of infor- 
mation about neutron structure function F 2 n (x). The main obstacle in the quantitative 
evaluation of F^x) in this case comes from nuclear binding which is neglected in many 
analyses on the grounds that the deuteron binding energy is very small. Of course, one 
can assume that it will be sufficient to consider just Fermi motion which becomes partic- 
ularly important at large Bjorken x. On the other hand the EMC effect [18] and nuclear 
shadowing [19] show us that even small binding can qualitatively change the observed 
nucleon structure. Its effect, as it is demonstrated by the analysis in Ref. [20], is clearly 
manifested in the entire x range, including the range of x < 0.1. Up to now there is no 
well established and unified explanation of nuclear effects in whole range of x. Numerous 
models used to reproduce dynamics of the effects are not consistent. All these facts render 
the procedure of F£(x) evaluation very ambiguous and model dependent. 
To elucidate the problem of model dependence, one needs an approach which is less 
dependent on the details of dynamical nature of the effects and provides more general and 
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unified picture. We use the approach based on the covariant Bethe-Salpeter formalism [13]. 
It yields a good description of the ratio of the nuclear to deuteron structure function which 
is shown to be universal for all nuclei [15, 20]. Within this approach the hadronic part of 
the nuclear deep inelastic amplitude W A V is expressed in terms of the off-mass-shell nucleon 

and antinucleon amplitudes, and W^ u , respectively, by the following expression: 



E / dk l (W^(k l ,q)f/ A (P,k t ) + W^(k l ,q)f/ A (P,k l )) , 



(1) 



where the indices [w denote the Lorentz components of the amplitude, index i counts 
nucleons inside the nucleus, P is the total momentum of the nucleus, ki is the relative 
momentum of the struck nucleon and q is the transferred momentum from the photon. 
The distribution function f N / A is expressed in terms of the n-nucleon Bethe-Salpeter 
vertex functions 

f/ A (P, ki) = Jdh... dki-idhi+i . . . dk n u%)S {n) {P, k t )u(^)7(P, £)S {n) (P, /C)r(P, /C), 
f/ A (P, ki) = fdh... dki-rfh+i ■ ■ ■ dk n v(ki)S {n) (P, ki)v(ki)T(P, K,)S {n) (P, K,)F(P, JC). (2) 



It gives AD momentum distribution of the nucleon and antinucleon inside a nucleus. Thus, 
according to Eq. (1), all nuclear effects should follow from the AD Fermi motion of the 
nucleon inside a nucleus [21]. The time component of the Fermi motion is exclusive feature 
of the relativistic approach. In the 3D limit this component results in the change of the 
nucleon structure [21]. Therefore, one can take explicit account of it by using a dynamical 
model for the interlinkage of the nucleon and nuclear structure. This problem is closely 
related to that of the off- mass-shell effects. Since the nucleon amplitude in Eq. (1) is 
off shell it cannot be connected with corresponding observable nucleon amplitude. This 
point makes Eq. (1) useless for practical applications. To solve this problem one can use 
analytic properties of the integrand in (1) and remove explicitly the Fermi motion along 
time axis by integrating over k . We do it assuming small relative momenta of bound 
nucleons. Preserving general form of Eq. (1) for the lightest nuclei one can transform it 
as follows: 



A-l 



W A (P,q)=E 



d? k a 



(ltv a Q 



(3) 



where a and a' denote the struck and spectator nuclear constituent, respectively. Depend- 
ing on the mass of the considered nucleus, a and a' assume different symbols: 



(a, a') 



(p,n),(n,p), if A = 2 (D) 
(p, pn), (n, pp), (p, D), (D, p), if A = 3 ( 3 He 
(n, pn), (n, pn), (n, D), (D, n), if A = 3 ( 3 H) 
(p, pnn), (n, ppn), (p, Dn), (n, Dp), 
(D,pn),(D,D),(p, 3 H),(n, 3 He), 



(4) 



if A = 4 ( 4 He) , 



where A A a , = Ma — E a — E a > is the separation energy of the corresponding nuclear 
fragment. Now the nuclear effects are interpreted by the conventional 3D Fermi motion 
of nuclear fragments and the derivative with respect to ko of its DIS on-shell amplitudes. 
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The distribution function $^ a /(k) 2 is defined by the projection of the Eq. (2) onto 3D 
space 

Ka^) 2 = {f /A (P,h)\ 2 (5) 

and it is closely related to the nuclear spectral function. The term containing the deriva- 
tive represents a dynamical modification of the bound nucleon structure observed in the 
3D projection of the relativistic bound state as nuclear shadowing and the EMC effect. 
The spectrum of bound states presented by Eq. (4) is defined by analytic properties of 
the distribution function f N ^ A and thereby ensures Pauli blocking of the forbidden states. 
Starting from Eq. (3) one can derive the F 2 D in the form [14]: 

D , N f d 3 k m 2 f^ N/ jE-h : M D -2E 



/ N /°(M D , fc) - F 2 N (x N )(£ - A; 3 W N/D (M D , fc) 



fc =-B-M D /2 



where E 1 is the on-mass-shell nucleon energy. We consider this expression as an integral 
equation with the unknown function F 2 n . Functions F 2 P and F 2 D are assumed to be known 
from experiments. Since the available data does not cover the entire kinematic range, 
the behaviour of these functions at boundaries (x — 0, x — 1) can not be defined from 
experiments. The latter makes definition of the boundary conditions for the solution of 
Eq. (6) model dependent. Basically, the boundary behavior of the ratio -F 2 n /F p can be 
found from an iteration procedure. Experience of using such procedure shows that the 
final result strongly depends on the zero iteration [12], which is hard to define on the 
basis of physical conditions. Another way of solving this problem is to fit the right hand 
side of Eq. (6) to experimental data on F 2 D (x). The boundary conditions are introduced 
explicitly as asymptotics of the nucleon structure functions at x = and x — 1 and their 
influence on F^i^x) at medium x is studied. 
As an initial condition we use the following anzats for F^x): 

F*(x) = R n ^(x)F^(x) (7) 

with the extended SMC fit of and the function R n / p (x) in the form 

R n ^(x) = a l {l-x) ai +a 2 x a2 +b 1 x^{l-xf 2 {l + c l x 11 ) . (8) 

The parameters a±, a 2 , «i, &i are introduced in order to satisfy the asymptotics of the nu- 
cleon structure functions: a\ — 1 (corresponds to the limit -F 2 P (0) = F 2 n (0)); the parameter 
a 2 should be fixed according to /im :E ^ 1 F 2 n (a;)/F 2 p (x). We use three values for this limit in 
order to study possibility to extract it from the experimental data. The simplest quark 
model with SU(6) symmetry gives i?°= p = 2/3. On the other hand, the kinematical limit 
x = 1 corresponds to the elastic scattering off the nucleon. Therefore, the ratio of the 
nucleon structure functions becomes equivalent to the ratio of the elastic cross sections - 
Rl-i = a D /a p .It gives another value for the ratio in this limit, i? 1 ^ = 0.47. 

x — i clastic' clastic J ^ — ^ 

The minimal value of the ratio has been provided by the model with SU(6) symmetry 
breaking with scalar diquark dominance — Rx=i = 1/4. Thus one has a wide range of 
possible values for the structure functions ratio at x — 1. In the next section we show 
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Figure 1: The structure function of the deuteron measured in the SLAC (filled squares) 
and NMC (empty circles) experiments, shown as a function of x for bins of fixed Q 2 (left 
upper panel — on a log scale, and right upper panel — on a linear scale), and as a function 
of Q 2 for bins of fixed x (two lower panels). It is approximated with Eq. (6) in the range 
10~ 3 < x < 0.6 with the constraints as explained in the text. 
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how measurements of FJ? / F% are sensitive to these assumptions. The last constraint for 
the parameters can be obtained from the asymptotics of the proton and neutron structure 
functions by assuming validity of the quark counting rule in the range x — > 1. This results 
in similar asymptotic behavior for F^(x) and F£(x): 

limF 2 p (x) ~ CJl -x) 3 , (9) 
limF 2 n (:r) ~ C„(l -x) 3 . (10) 

x— >1 

Accordingly, the derivative of i? n / p at x — 1 is zero because 

lim R n / p (x) = -± = Const. (11) 

X^l Cp 

This gives following constraints on the parameters of Eq. (8): oc\ — 1, fa — 1, &i = 
(01202 — 1)/(1 + ci). All other parameters are considered as free and used to fit Eq. (6) to 
the deuteron data in the range 10~ 3 < x < 0.6. 



3 Discussion of results 

The procedure described in the previous section is used to approximate the SLAC and 
NMC data on F 2 D (x, Q 2 ) in the range 0.6 GeV 2 < Q 2 < 65 GeV 2 and 10~ 3 < x < 0.6. This 
is done by fixing three different limits of F^jF^ at x — 1 and by varying four parameters 
in Eq. (8), namely 0-2, fa, c\ and 71. In the considered kinematic range all three limits 
provide equally good approximation, which means that the solution for F^[x) j ' F%{x) 
converges in the range x < 0.6 to the virtually unique function described by Eq. (8). 
The result of the fit, which corresponds to the case F 2 n (l)/F 2 p (l) = 0.47 is displayed in 
Fig. 1. The three alternative solutions are shown in Fig. 2 with three lines, which virtually 
coincide if x < 0.4. Outside this range, the obtained ratio F^x) / F^{x) is represented by 
three lines, which exactly reproduce the constraints imposed at x — 1. The convergence 
to three different solutions is apparently provided by two constraints: available data on 
F 2 D and by three considered boundary conditions. Therefore, the solution is always model 
independent in the range of x, in which the constraint from F 2 D measurement is stronger 
than the one from boundary conditions. 

The solutions, which we find for F£(x) are compared in Fig. 2 with the results of NMC 
experiment [22] and SLAC analysis [23]. The NMC results are obtained assuming a naive 
approach F 2 n = 2F 2 D — F 2 P that is equivalent to the assumption F 2 D = F 2 N . The points 
corresponding to SLAC results were transformed by us from F®(x) / F^x) to F^x) / F^x) 
in the same naive way as it was done by NMC, namely by neglecting nuclear effects in 
the deuteron. In the range x < 0.6 all three lines are in good agreement with the results 
of NMC and SLAC. We explain the success of the naive approach in the range x < 0.6 by 
the effect of cancellation of the modifications of the free nucleon structure function due to 
the nuclear (deuteron) binding and nucleon Fermi motion in this very kinematic range. Of 
course, the cancellation is not complete, which is seen from small (2-4 %) overestimation 
of the data at medium values of x. The difference between the naive approximation and 
the exact result becomes dramatic at 0.62 < x < 1, which is a well known problem. This 
is well illustrated in Fig. 2 by inconsistency between data points at x < 0.8 and x > 0.8. 
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Figure 2: The ratio of the neutron to proton structure functions evaluated in the frame- 
work of the presented approximation by setting three different values for the parameter 
R^i- (1) — 2/3 (short-dashed line), (2) — 0.47 (solid line) and (3) — 1/4 (long-dashed 
line). Points with error bars depict the results of the NMC experiment [22] (solid squares) 
and SLAC analysis [23] (solid triangles) obtained on the basis of the naive approach. Ad- 
ditionally, the small x range (x < 0.1) is displayed on a log scale in the inset. 

To restore the consistency, complete accounting of both the nuclear binding and Fermi 
motion has to be done in evaluating R n / p (x). 

In order to understand theoretical uncertainties in description of the data one should 
investigate how the extracted F^x) can modify theoretical evaluation of the deuteron to 
proton ratio. The three different constraints on R^=\ considered in this paper change the 
approximation of F^(x) in the range of high x which is better seen in the ratio F 2 D /F 2 P 
shown in Fig. 3, in which three alternative calculations are compared with data from 
Ref. [22] and the results of the analysis performed in Ref. [23]. We conclude, that the 
present status both of the data and theory does not allow to constrain the value of the 
function R n lv{x) at x = 1: both 1/4 and 0.47 are in agreement with experiments, while 
the value 2/3 can be regarded as less preferable. Further understanding of the relation 
between F^x) and F^{x) can not be achieved without improvement of experimental data 
on F® (x) / F£ (x) in the range above x = 0.6. The latter is also valid for points from 
Ref. [23] even if they have apparently smaller error bars. Contrary to the NMC data 
which have similar x dependence in different Q 2 intervals, the data from SLAC [24] used 
for the analysis [23] scatter considerably, which is not consistent with the parton picture 
suggesting the valence quarks dominance in the high x region. Our understanding of the 
problem is that the NMC data comes from the single experiment, whereas SLAC results 
are based on eight experiments on deep inelastic e — p and e — d scattering. 
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Figure 3: The result of modification of the ratio of the deuteron to proton structure 
functions evaluated for three values of the parameter R^.=i as in Fig. (2) compared to the 
data of Ref. [22] (solid squares) and results of the analysis performed in Ref. [23] (solid 
triangles). Additionally, the small x range is displayed on a log scale in the inset. 



We would like to note that the analysis performed in this paper by employing Eq. (3) is 
consistent with the data available for the ratio of 4 He and deuteron structure functions in 
the range 10 -3 < x < 0.8 [13, 15]. Good agreement with experiment in two different cases 
observed in a wide x range proves that Eq. (3) offers a general approach for accounting 
of the nuclear effects in the whole region of x. 

Naturally, F 2 D (x) can be further modified in the high x region if one assumes the presence 
of non-nucleonic degrees of freedom in a nucleus, which are not implied in the nucleon 
structure. Moreover, if there exist dibaryon states, like 6q states, which can not be 
excluded on theoretical grounds, F 2 T) (x) can significantly change at x ~ 1. Our calculations 
therefore provide reference lines for the search of such effects in forthcoming measurements 
proposed for the upgraded CEBAF facility [25]. 



4 Conclusions 

We have proposed theoretically justified and fully consistent procedure for extracting 
the neutron structure function F 2 n (x) in the kinematic range 10~ 3 < x < 1 under three 
different assumptions on F 2 /F 2 at x — 1. The procedure involves a numerical fit of the 
expression (6) to the deuteron data. The performed analysis indicates that the increase 
in experimental accuracy in measurements of F 2 D (x) / F^x) in the range 0.6 < x < 0.8 
by factor of two will be sufficient for verification of models suggested for the evaluation 
of the d/u ratio at x = 1. The developed procedure allows one to avoid appreciable 
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theoretical ambiguities which are present in other analyses largely due to simplifications 
in the treatment of Fermi motion. This concerns rather wide interval of x and not only 
high x range as it is commonly believed. The procedure proved to be a robust one 
because it relies on a good approximation of F^(x) which is not sensitive to different 
high x limits of the neutron structure function. This also means that F®(x) measured 
by already completed DIS experiments (x < 0.9) can be described without introducing 
nonbaryonic degrees of freedom. The interval which remains unmeasured can in principle 
accommodate dibaryon states or some other exotica. More data in the high x region is 
required in order to obtain model independent information on hadronic structure of the 
deuteron and to find which physics is responsible for the d/u ratio. 
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6 Appendix 

Here we present for completeness the parametrization for F$(x, Q 2 ) as suggested in 
Ref. [17] and slightly modified in the present analysis: 

// (Ql)\ B{X) I £ K nX n \ 

if(*,Q 2 ) = ^(i-*) A2 £ c n (i-x)^ lp$>) 1 + n=1 "« (A.l) 

n=1..5 \Mf)/ V Q J 

where 

B(x) — pi + p 2 x H — — . 

Fit parameters as obtained in Ref. [17] are presented in Table 1. The parametrization 
is restricted to the kinematic region 3.5 • 10~ 5 < x < 0.85. In order to extend it to the 



% 


\ 


Pi 




a 


1 


-0.2499713175097 


0.1141083888210 


-1.451744104784 


0.2289630236346 


2 


2.396344728724 


-2.235597858569 


8.474547402342 


0.08498360257578 


3 




0.03115195484229 


-34.37914208393 


3.860797992943 


4 




0.02135222381130 


45.88805973036 


-7.414275585348 


5 








3.434223579597 



Table 1: Values of the parameters for F 2 P given in Eq. (A.l). 



region of high x and satisfy the (1 — x) 3 behavior at x — > 1 we modify the parameter A 2 
as follows: 

A 2 -> A 2 (x) = A 2 + (3 - A 2 )x 15 . 
9 



This correction does not affect the values of F 2 at x < 0.6 and affords an approximation 
of the proton data in a much wider kinematic region. 

The neutron structure function is defined by Eqs. (7) and (8) in the main text. Taking 
into account the constraints on the parameters we arrive at the following expression for 
i? n / p : 

R*/p(x) = (l-x) + a 2 x a2 + a2Q2 ~ 1 x fh (l - x)(l + c x x^) . (A.2) 

1 + C\ 

The fit parameters are presented in Table 2. 





a 2 = 2/3 


a 2 = 0.47 


a 2 = 1/4 




3.13971 


2.2262 


1.15416 


Pi 


2.2129 


1.61188 


0.88126 


Cl 


-1.01176 


-1.00692 


0.86217 


7i 


0.01901 


0.08483 


5.65744 



Table 2: Parameters of Eq. (A.2) which is the constrained form of Eq. (8) connecting the 
proton and neutron structure functions. 
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